EXISTENCE AND NON-EXISTENCE OF GLOBAL SOLUTIONS FOR 
A DISCRETE SEMILINEAR HEAT EQUATION 
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Abstract. Existence of global solutions to initial value problems for a discrete analogue of a 
d-dimcnsional semilinear heat equation is investigated. We prove that a parameter a in the partial 
difference equation plays exactly the same role as the parameter of nonlinearity does in the semilinear 
heat equation. That is, we prove non-existence of a non-trivial global solution for < a < 2/d, and, 
for a > 2/d, existence of non- trivial global solutions for sufficiently small initial data. 
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1. Introduction. The blowing up of solutions to the semilinear heat equation 
has been analysed extensively since the pioneering work by Fujitap] |2 [3]. Fujita 
studied the initial value problem of the equation: 




A/ + / 1+ « (a>0) 
) = a(x) > (a(x) ^ 0) 



where / := f(t, x) (i>0,3?€ R d ) and A is the d-dimensional Laplacian A := ^ -J^-. 

fc=i Xk 

When the initial value a(x) is continuous and unifomly bounded, there is a smooth 
solution for t > and whenever the solution is bounded, the solution is prolongable. 
Moreover, since a(x) > 0, the solution satisfies f(t,x) > 0. A feature of (11.11) is that 
the solution is not necessarily bounded for all t > 0. This fact is easily understood 
if one considers spatially uniform initial condition, a(x) = a G M + . In this case, 
f(t,x) = f(t) and (jl.lj) becomes an ordinary differential equation, 

f 4f _ f l+a 
(1.2) Idt J 

[/(0)=a>0 

The solution to this differential equation is 

f{t) = (a-i a -«-t)V«' 

and we see that it diverges as t — > a _1 a _Q — 0. In general, if there exists a finite time 
T e R+ and if the solution of CO]) in {t, x) e [0, T) x R d satisfies 

limsup ||/(t, = oo, 

where 

||/(t,-)IU- := sup \f(t,x)\, 
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then we say that the solution of (jl.ll) blows up at time T, and therefore, that it is not 
a global solution (in time). In 1966, FujitaQ] proved the following theorem 
Theorem 1.1. 

(1) If < a < 2/d, any solution to is not a global solution in time. 

(2) If 2/d < a, then global solutions to flU!} do exist for sufficiently small and 
smooth initial functions a(x). 

A remarkable point in this theorem is that the parameter a affects the existence 
of the global solution of (jTTTJ) . The critical value 2/d in this theorem is called the 
Fujita exponent. The case a = 2/d was studied by Hayakawa for d = 1,2, and by 
Kobayashi, Sirao, Tanaka[5], and Weissler[6] for general d. They proved, 

Theorem 1.2. For the case a — 2/d, there exists no global solution to . 

In numerical computation of (jl.ip , one has to discretize it and consider a partial 
difference equation. A naive discretization would be to replace the i-differential with 
a forward difference and the Laplacian with a central difference such that turns 
into 

/•r+1 _ Pf d r r _ Q fT _i_ fT 

■In J_n_ _ \ * Jn+Sk Jn' Jn—e k , l rr\l+a 

k=l ^ 

where /(r, n)(=: /I) : Z>o x 7L d — > R, for positive constants S and £Cand where 
e*fc 6 Z d is the unit vector whose fcth component is 1 and whose other components are 
0. Putting A := (5/£ 2 , we obtain 

(1.3) fZ +1 = 2d\M(fZ) + (1 - 2dX)fl + S(fZ) 1+a (a > 0). 

Here 

1 d 

(1-4) M(Vn) := Yd E + V n-^)- 

k=l 

For a spatially uniform initial condition, (|1.3[) becomes an ordinary difference equation 

r +1 = r + s(D 1+a - 

The above equation is a discretization of (|1.2I) . but the features of its solutions are 
quite different. In fact, f T will never blow up at finite time steps. Hence, (|1.3p does 
not preserve the global nature of the original semilinear heat equation (ll.l[) . 

In this article, we propose and investigate a discrete analogue of (jl.ll) which does 
keep the important characteristic of existence and non-existence of the global solutions 
in time. In section 2, we present a partial difference equation with a parameter a 
whose continuous limit equals to (11.11) , and state the main theorem which shows that 
this difference equation has exactly the same properties as (II. ip with respect to a. 
This theorem is proved in section 3 for the case < a < 2/d. In section 4 for 2/d < a 
and in section 5 for a = 2/d. 

2. Discretization of the semilinear heat equation. We consider the follow- 
ing initial value problem for the partial difference equation 

|V+i gfl (r e Z> , ne Z d ) 

(2.i) < {i-{gm i/a - 

W = a(n)>0 (a(n)^0), 
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where a > is a parameter and gZ is defined by means of M (jl.4l) as 

By scaling / with a positive parameter 5 as 

FZ := (aS)-^fn, 
GZ := M{FZ), 

we have 

GZ 



(2.2) = 



{1 - a5(GZ) a } 1 / c 



If there exists a smooth function F(i, x) (t S R>o, £ € ]R d ) that satisfies F(tS, £ri) 
FZ with £ := V2dS, we find 

F(i + 5, x) = G(t, x){l + S(G(t, x)) a } + 0{5 2 ), 

with 

1 d 



2d 

k=l 



or 



F(t + 6,x)- F(t, x ) _ A F(t, x + £e fc ) - 2F(t, x) + f - je k ) 



2 

l+a 



+ (F(t,x)y +a + 0(6). 
Taking the limit S — > +0, we obtain the semilinear heat equation (jl.ip 

BF 

= AF + F 1+a . 

at 

Thus (|2.1[) can be regarded as a discrete analogue of (II. ip . 

Because of the term (1 — (#5)) , if — > 1 — 0, then /I — > +oo, and (|2.1[) cannot 
be defined when gZ > 1 for generic a. This behaviour may be regarded as an analogue 
of the blow up of solutions for the semilinear heat equations. Thus we define a global 
solution of (I2.1[) as follows. 

Definition 2.1. When a solution fZ of (|2.1|) is non-negative and uniquely 
determined for all r € Z > o and n € Z d , i.e. gZ < 1 for all r G Z 
i/ien we saj/ i/iat the solution fZ is a global solution (in time) of (|2.1|) . 

The advantage of using (|2.1[) instead of (( 1 . 3|) is apparant from the following propo- 
sition. 

Proposition 2.2. No spatially uniform function can be a global solution of 
CO). 

Proof. When the solution is spatially uniform, then (|2.1[) takes the form of an 
ordinary difference equation 

CT+l f T 



r +i = 



{i-(/ r )n i/Q 

{(f T ) a - l}- 1/a (for / + 0). 
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Hence, if a non-trivial global solution exists, it satisfies < f T < 1 for all r e Z >0 
and 

tr +1 r a = (rr a - 1. 

However, the above equation is easily solved and we get 

(rr a = (fr a - r 

Since (f°)~ a is a positive constant, f T cannot be defined for r > (f°)~ a , which is a 
contradiction. □ 

Furthermore, (I2.f [) inherits quite similar properties to those of The follow- 

ing theorem is the main result in this article. 
Theorem 2.3. 

(1) For < a < 2/d, there is no global solution to §2. 

(2) For a = 2/d, there is no global solution to 

(3) Let ||/° ||i be the I 1 norm of the initial function, i.e., \\f°\\i '■— X)« I /si- For 
2/d < a, if ||/° ||i is sufficiently small, then global solutions to $2.1$ exist. 

Remark. • In the limit a — > +0 in (12.21). we have 




= a{n) > (a(n) 0) 

Although || F T ||i diverges as r — > +oo as far as 5 > 0, the solution of this difference 
equation is a global solution in time. On the other hand, when a = 0, the partial 
differential equation (jl.ip becomes linear, and its solution is also a global solution. 
• To consider negative solutions or oscillatory solutions, we have only to use a slightly 
modified partial difference equation 

rr + l 9ti 



{i-(\g T n\) a } 1/a ' 

3. Proof of theorem 12.31 for < a < 2/d. The idea of the proof in this 
and the following sections is similar to that adopted by Meier in the case of partial 
differential equations [7]. We construct a subsolution and a supersolution of /I and 
prove the existence and non-existence of the global solutions. 

We denote by U~ the solution to the initial value problem of the linear partial 
difference equation 



(3.1) 




Mm) 



Using [/I, we define 

(3-2) ^ := ^ [/!_ .,/«,, 



and 

(3-3) fl:= 
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provided that 1 - T(hZ) a > 0. 

Remark. • Since the support of U T is finite, the summation in the definition of 
hZ is over a finite number of lattice points n. 

• Due to the definition of UZ, it holds that 1 1 7^.^~ 1 1 x = ||/°|| r 

• The function hZ satisfies the linear partial difference equation 



(3.4) 



= M ^l) 
h°n = f°n 



• Since hZ > 0, fl > if it exists. 

— — n 

• When hZ > and 1 - r{hZ) a > 0, 

(3.5) (flY a = (hZ 



• The function f T _ does not always exist for all r € Z > q. When 1 — T(hZ) a < 0, f T _ 
does not exist. 

To understand the meaning of fZ and in order to give the proof of theorem l2.3l f 1). 
we need several lemmas and propositions. First we consider the following elementary 
inequality. 

Lemma 3.1. Let Xk (k = 1,2, N) be non-negative real numbers, a > 0, and 

N 



Xk > 0. It then holds that 



k=l 



k=l 

N x 1+a 

k=l 



> 1. 



Proof. If we put 



N I N \ 

[x x , ■ ■ ■ , x N ) := N a £ (^) 1+Q -($>*) 



l + Q 



lemma [3Tl is equivalent to the inequality 4>{x\, ■■ ■ ,xn) > 0. Without loss of gener- 
ality, we can assume x\ > X2 > ■ ■ • > xn > 0. From the partial differentiation 

d 

"xi,--- ,x k ,x k+ i,x k +i, ■ ■ ■ ,x k+ i) (k = 1,2, ... , JV - 1) 



N-k 

k 



we find 



= (JV-fc)(l + a)^ (JVa: fc+ i) Q - ( E^ + ^-fcK+i ] }> <0, 

^=1 



(Si, • • • ,Xk,Xk+l,Xk+l, ■ ■ ■ ,Xk+l) > <f>(xi, ■ ■ ■ ,Xk-l,X k ,X k , ■ ■ ■ ,x k ). 

V v ' " * ' 

N-k N-k+1 



G 
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Hence we have, recursively, 

<f>(xi, ■ ■ ■ ,XN-!,Xn) > <f)(x-L,- ■ ■ , X N -2,X N -i,XN-l) 

> (()(X1, ■ ■ ■ ,X N - 3 ,X N -2,X N -2,XN-2) 

>•■•> 4>{xi, ■■ ■ ,xi) = 0. 

□ 

Lemma 3.2. Suppose that fZ, and /^ ±e - (k = 1,2, ... ,d) exist and that M(fZ) > 
0. 7/{M(/T)}~ Q — 1 > 7 then exists and satisfies 

(3.6) (/r 1 )-" > {M{QV a 1. 

Remark • The condition M(fZ) > implies that at least one /^ ± - is positive and 
that hl +1 = M{hV) > 0. 

• When M{h T R ) = 0, then hl +1 = and fl +1 = 0. 

Proof. We assume that there are m (1 < m < 2d) nonzero values among 
{hft±e k }k=i: an d denote them by hi,h,2,---, h m . Accordingly we put 

fi ■= tt = (K a ~ T )~ 1/a (i = 1,2,..., m). 

(l-r(/, 4 )") 1/Q 

By definition, 

X (r) :={MK)r-r-{M(/l)r 

= { Yd 2 + J } - T - | ^ £ + ) } 



Regarding r as a continuous variable, 



dxw _ . 1 + (2dr e::i{(^)- q -^}-( q+i )/ q 



rfr [E£i{(^)- a -r}-V«] aH 



> -1 + (m) c 



K"=i{('if)- a -t 1/o ]' 



Appling lemma O for N = m, x t = {{h t y a - t}" 1 /", we find > o. Thus 

GIT 

x( T ) > x(0) = 0, and we obtain 

[M{hl)}~ a -r>{M{fl)y a . 
By the assumption {M[fV)}~ a - 1 > and (|3~3|) we have 

- (r + 1) = {M(^)}- a - r - 1 > {M(/D}-° - 1 > 0. 
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Thus we find 

i-(T+i)(/£ +i r>o. 

Therefore exists and due to (|3.5[) it satisfies 

(r H +1 r a > {M(f r )r a - 1. 

□ 

Proposition 3.3. If the solution fZ of the initial value problem (|2.1[) exists at 
t and for all n € Z d , then f_ exists at the same r and for all n € 1i d , and satisfies 

(3.7) fl < f T n . 



Proof. The proof goes by induction for r. Since /° = /2 by definition, fZ 
exists and satisfies (|3.7j) for r = and for all n. Suppose that fZ exists and satisfies 
(ET7| for all n and for all r < s. If exists, either = or > 0. When 
/2 +1 = 0, it implies 

fl +1 = < 



M(fn) 


i = 








J n±ek 


= 


(k = 


1,2,.. 


,,d) 


f 

—n±e k 


= 


(k = 


1,2,. 


..,d) 


n n±e k 


= 


(k = 


1,2,. 


..,d) 



M(h s H ) = 
= 

n 

^ ft +1 = o. 
—n 

Hence fl +1 exists and satisfies ff, +1 < f i + . 

When ff +1 > 0, if M{ft) = 0, then M{hi) = and = 0. Otherwise 

{M(Ln)}~ a ~ 1 ^ W^}"" - 1 = (fn +1 T a > 0. 
Then, from lemma l3~2l the existence of follows. Moreover, from (|3.6[) . 

(/f T" > {M(/^)}- Q - 1 



r-s+l ^ JS+l 

— n — ■'ft 



Thus we find 



From the induction hypothesis, (|3.7p holds for arbitrary r 6 Z+. □ 

The following asymptotic evaluation of UZ is well known from the analysis of 

the transition probability of d-dimensional simple random walk. See for example 

F. Spitzer[5]. 

PROPOSITION 3.4. 

/ j \ d/2 

(3.8) UZ~2\J-\ T - d /2 {T ^ + ^y 
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Here A T ~ B T (r -> +00) means lim T _j. +00 (A r /_B T ) = 1. 
Now, we are ready to give the proof of theorem 12.31 (1). 
Proof of theorem\EE (!)■ 

Suppose that the global solution fZ exists. From proposition 13.41 is evaluated 

as 



hZ = X^ IJl f 2 

n / j n—n' J n 



(3.9) ~vr^ 2 {i) $ {T ^ +oo) - 

n 

Putting C :=2(d/4%) d/2 \\f\\ 1 , 
(3-10) fl 



-H _ Qa T l-da/2\l/ a ' 



For < a < 2/d, the exponent of r in the above equation satisfies 1 — da/ 2 > 0. 
However, since C is a positive constant independent of t, 1 — c a T 1 ~ da t 2 < for 
sufficiently large r and the evaluation (|3.10|) does not make sense. In other words, 
fZ cannot exist for sufficiently large r, which contradicts proposition 13.31 and the 
statement (1) of theorem 12.31 is therefore proved. □ 



4. Proof of the theorem 12.31 for 2/d < a. First, we define a supersolution 
of JO 



(4-1) II:-- t ''■ 

' 1 - E M a 

k=0 



1/a ' 



where m T is defined in terms of (|3.2[) as 

(4.2) m T := sup h T n . 



Of course fZ is well defined only when 1 — / jQrijfe)" > 0. 

fe=0 

Proposition 4.1. Mien /I exists at r and for all n £ Z d , /I exists and satisfies 
(4-3) i£ > fh 



Proof. We give the proof by induction on r. When r = 0, by definition of the 
initial value problem, /2 exists and (I4.3[) holds because 
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Suppose that the statement is true up to r = s and exists. When = 0, 

J n, 











M{h s n ) = 








h n±e k = 


(k = 


1,2,. 


..,d) 


fn±s k = 


(k = 


1,2,.. 


..,d) 


fn±s k = 


(k = 


1,2,.. 


..,d) 


.9^ = 








4 2+1 = o. 









Hence (jO|l holds. 

When fS +1 > 0, if q- = 0, then fi +1 — and the statement is true. Otherwise 

i - E (™fe) Q i - E / \ « 

n < (f s + 1 \~ a = fc =° = fc =° ( m s+i \ 

1 - E (m k ) a 



< (g s n)- a i- 

From (EH), - 1 = and we find < i.e. < 

Thus, from the induction hypothesis, the statement is true for any non- negative 
integer r. □ 

Now we prove the statement (3) of theorem 12.31 
Proof of theorem\2J\ (3). 

From (|3.9p . we obtain the asymptotic behaviour of m T as 
m T ~ — — [t — > +oo j. 



Hence with a fixed ro € Z>o, we can evaluate /I as 

Zi- 
yi 12 ( T 

1 - E m fc - E 7^ 



oo 1 

Since da/2 > 1, we have ^ TdaH < Noticing the fact \\h T \\ l = \\f°\\ v 



k da/ 

fc=TQ + l 



the term E m fc an d the constant C can be as small as possible by choosing a small 

fc=0 

value of H/ )^- Thus fZ exists at any time step r and any lattice point ft, if H/ ^ is 
sufficiently small. Then, from proposition ^. 1[ we find that fZ is a global solution of 
(I2~T1) . □ 
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5. Proof of theorem 12. 31 for a — 2/d. In this section, we prove the statement 
(2) of 12.31 The idea of the proof is similar to that adopted by Weisslerj^. 
We define the discrete Green function G and the function / as follows. 



(5.1) GZ 



(t = 0) 
UZ- 1 (r>l) ' 



(5-2) fZ r- 



(r = 0) 

s=0 n' 



where 

(5.3) H(g) :-- 



(1 - g a y/ a 



Since 



fr+l _ T _ 9n t TS( n T \ 



the initial value problem (12.11) can be rewritten as 
(5.4) 



(T — M)fZ = H(gZ) 
.§=a(n)>0 (a(n)#0), 



where T denotes the time shift: Ty T := y T+1 . 
For r > 1, 

(f - M)GZ = GZ+ 1 - M(G T R ) = UZ - M(UZ^ 

and for r = 0, 

(f - M)G% = G\ - M{G%) = C/2 - ^, . 

Thus we find 

(5.5) (f-M)GZ. = 5 Tt0 5 flt0 . 

We also find 

(5-6) (T-M)fZ = H(gZ). 

Because 

(f - M)fl = }Z +1 - M(fZ) 



E E G r+-* S H(9k>) + E E 

s=0 n' s=0 n' 

T-l 

J2 G\^H(gZ,) + E E ( f - M){ GT n-% H ^)} 



s=0 

T-l 



E u n-rvH(gZ,) + E E {*r-.,o**-ff 

n' s— n' 

H{gZ). 
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Using the above results, we can prove the following proposition. 
Proposition 5.1. When fZ is a global solution of (|2.ip 

T 

(5.7) « = ^+EEw^ 



Here hZ is defined in 

T 

Proof. First we note that fZ = \^ UZ~i,H (giV 1 ), and the right hand side of 

s=l n' 

(T5T7T) = hZ + fZ. By definition h% = /9 and f% = 0. From (gUJ and 



[(T-M)(hZ + fZ)=H(gZ) 

Thus, from (|5 .4[) . hZ + fZ satisfies the same initial value problem as that for fZ. Since 
the solution of the initial value problem (|5.4p f (|2.ip ) is unique, fZ = hZ + fZ. □ 

PROPOSITION 5.2. Suppose that a = 2/d and that fZ is a global solution of (|gjp . 
Then there exists a constant Co and 

ll/li<C . 

/or a/Z r G Z + . 

Proof. From proposition 13 .31 and the assumption of the proposition, fZ, also exists 
for all (r, n) and 

(5.8) 1 - {T d/2 h T n ) a > 0. 

On the other hand, from proposition 13. 4[ 

d/2 / d \ 

4tt ) fn = 2 l 47 J 



Then from (Oil . 



Therefore, 



»i£)*Vi.} <l 



d/2 



Because /I is a global solution, if we take /I at any time step r as the initial value 
of (|2.1j) , the solution is also a global solution. Therefore we obtain 

ll/l|l<: 2UJ ' 

□ 

Lemma 5.3. When fZ is a global solution, 
(5.9) ^ < /I. 
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Proof. We prove by induction for r. 

h°- = V C/2 „. f° = f 2 

ft' 

and (|5.9p holds for r = 0. Suppose that the statement is true for r = s. then 

p+i = M(/4) 

" / . \l/a 

(l-M(/£)« 



(l - M(/$)° 
> M{K) = hi +1 

Here we used the inequality 

< x < y < 1 =^> ^ , . < 



l/a 



{l-x a ) 1/a ~ {l-y a ) 1/c 



Thus, by the induction hypothesis, (|5.9[) holds for any r £ Z > o- D 

Now we prove the statement(2) of theorem 12.31 . 
Proof of theorem\EM (2). 

Suppose that fZ is a global solution of (|2.ip . From proposition 15.11 

iiriii^EE^-S'W 1 ) 
= EE^r 1 )- 

Here we applied yi~ UZ = 1. By an elementary computation we can easily show 

(5.10) H(g)>-g l+a (0 < g < 1). 

a 

Hence 

(5.H) \\r\\i> EE -for 1 ) 1 -**- 

s=l n 

From (1531) . 

it' 

There is at least one no £ Z d such that f% Q > 0, and we put c := /2 Q > 0. Then 

^>c[/2_ llo , 
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and 

ii/i^EE^) 



(5-i2) ^ cl+Q -EE(^) 



Moreover, from proposition 13.41 we obtain the asymptotic behaviour 

(5.13) U *~ 2 (i) dl2 ^ d ' 2 ^( J €) (T ^ + ^ 

Let 8 and £ be 5 := — and £ := \/2dS respectively. From ()5.13p . we obtain the 
r 

following evaluation for a = 2/d and r — > +oo (5 — > +0, £ — > +0). 

? ™ w ' - E 2 ' +2,J (sf) "° »>' {-f (' + 1) } +~> 

71 ft 

2 i + mJ V exp / -« fl + AU 

4^r (87rr(5) d / 2 ^ P \ 8tS \d + d 2 Jj 

_d2^/i _2\- d/2 i 

oo r 

Since X) V s diverges, Yl (Ufi) 1+2 ^ d can take an arbitrarily large value. From 

s— 1 s— 1 n 

(|5.12p . ||/ r |ji also becomes arbitrarily large, which contradicts with proposition 15.21 
Therefore, when a = 2/d, there exists no global solution and we have completed the 
proof of theorem 12.31 □ 
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